A NON-ITERATIVE ALGEBRAIC RECONSTRUCTION 
TECHNIQUE FOR TOMOSYNTHESIS 

TECHNICAL FIELD 
[0001] The present invention relates generally to digital imaging, and 
more particularly to reconstructing a three-dimensional (3D) image from 
projections acquired using a tomosynthesis device. 

BACKGROUND OF THE INVENTION 
[0002] Tomosynthesis is an imaging modality where typically 
projection radiographs from only a few angles within a relatively narrow 
angular range are acquired. From these projection radiographs, a 3D 
representation of the structure of the imaged object can be reconstructed. 
However, since only "incomplete" information exists (i.e., not densely spaced 
projections over the full angular range), the reconstruction problem is difficult 
to solve. 

[0003] The main problems that need to be addressed by a reconstruction 
algorithm in order to yield satisfactory image quality are (i) efficient separation 
of overlying tissue, (ii) enhancement of the contrast, particularly of small 
structures, and (iii) minimization of artifacts. In some advanced reconstruction 
algorithms, some form of re-projection consistency constraint is utilized 
(directly or indirectly) to obtain high-quality reconstructions. These algorithms 
include, for example, additive algebraic reconstruction techniques (ART), 
matrix inversion tomosynthesis (MITS), Fourier based reconstruction, and 
volumetric order statistics-based reconstruction. However, these algorithms are 
generally computationally intensive and/or not very flexible to use. For 
example, strategies for artifact management may be hard to integrate in these 
known algorithms of the prior art. 

[0004] A need exists for an improved, computationally efficient and 
flexible method for reconstructing a 3D image of an object from a plurality of 
projection radiographs. 



SUMMARY OF THE INVENTION 
[0005] The present invention comprises a method of computing a 
reconstruction of the imaged volume that is equivalent to (in one embodiment), 
or approximates (in another embodiment) a reconstruction obtained by additive 
algebraic reconstruction techniques (ART). Generally, ART requires several 
iterations before a good image quality is obtained. In contrast, with the present 
invention, no backprojection/re-projection iteration of the dataset is required. 
The iterative update is replaced by a preprocessing step (filtering and linear 
combination), followed by a single backprojection step, thereby creating 
considerable computational savings. Furthermore, since the method does not 
involve a re-projection step, it is not required to reconstruct the full volume in 
the backprojection step, but instead the focus can be directed to the volume of 
interest, which may be a small region of the full imaged volume. In addition, 
the computational savings of the inventive method as compared to standard 
ART reconstruction allows the incorporation of more advanced processing 
schemes, while computing a reconstructed 3D representation of the imaged 
object in the time allotted, thus further improving the reconstructed image 
quality. For example, the present invention can be combined with additional 
artifact minimization (e.g., OSBP — order statistics-based backprojection) or 
contrast enhancing methods, while still maintaining a reasonable level of 
computational effort. 

[0006] In accordance with the present invention, projection images are 
acquired of a volume of interest (which contains, for example, portions of an 
object or person) at a plurality of projection angles. Filtered versions of all 
projection images are then prepared and combined in a suitable way to from 
new images and the 3D imaged volume of interest is then reconstructed from 
the new images. One aspect of the present invention covers the design of 
optimal filters for this reconstruction approach. 

[0007] Through the flexibility of the present method, a number of 
modifications can easily be derived, such as, for example, modifications that 
lead to reduced noise effects, or creating a filter for a "traditional" filtered 
backprojection type reconstruction. 
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[0008] The present invention is particularly useful, for example (but not 
by way of limitation) for producing reconstruction images of high quality in the 
context of tomosynthesis for chest and lung imaging. Numerous other uses of 
the present invention are possible as will be understood by persons of ordinary 
skill in the art. 

[0009] The above and other objects, features, and advantages of the 
present invention are readily apparent from the following detailed description of 
the preferred embodiments for carrying out the invention when taken in 
accordance with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] FIGURE 1 shows the basic configuration of an exemplary 
tomosynthesis system. 

[0011] FIGURE 2 illustrates the shift in the image of the center slice of 
the imaged volume with respect to different focal spot positions of the x-ray 
source. 

[0012] FIGURE 3 illustrates the situation where the 
backprojection/reprojection can be approximated by a convolution with a 
boxcar filter. 

[0013] FIGURE 4 illustrates examples of filters for a reconstruction 
from five projections. 

[0014] FIGURES 5 illustrates filters derived from the filters depicted 
from Figure 4. 

[0015] FIGURE 6 illustrates the derivation of an exact filter that can be 
used as a building block to determine the filters used in the reconstruction. 

[0016] FIGURE 7 is a flow diagram of another embodiment of the 
invention. 

DESCRIPTION OF PREFERRED EMBODIMENTS 
[0017] Digital tomosynthesis is a three dimensional (3D) x-ray imaging 
technique, where typically only a few projection radiographs are acquired for 
varying positions of the x-ray tube with respect to the imaged object and the x- 
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ray detector. Figure 1 depicts an exemplary tomosynthesis system 10 used in a 
wall stand system for chest/lung imaging. In particular, the detector 12 and the 
imaged object (i.e., the patient) 14 are stationary, while the x-ray tube 16 is 
movable along a linear trajectory 20 at a constant distance 18 from the detector 
plane, in order to acquire projections from different view angles. The detector 
consists in this example of 2048 x 2048 pixels, with a pixel size of 200 x 200 
jim, which translates into a detector size of about 41 by 41 cm. The tube- 
detector distance 18 (at the central tube position with respect to the detector) is 
about 180 cm, and the total range of motion of the tube is generally between 
about 31.5 and 131.0 cm. This corresponds to an angular range of ±5 to ±20 
degrees, where 0 degrees corresponds to the central position of the tube. 
Typically at least ten projection radiographs are acquired during a 
tomosynthesis acquisition, for tube positions covering the full angular range. 
[0018] Although the present invention will be discussed relative to the 
system 10 shown in Figure 1, it is to be understood that the present 
tomosynthesis reconstruction method is not constrained to this specific system, 
and can be used in other situations where the goal is to reconstruct 3D 
information about the imaged object/anatomy from relatively few projection 
radiographs. In particular, the detector specifications, as well as its position or 
orientation, may be different. Also, the detector may move during the 
acquisition, the x-ray tube may be moving along a non-linear trajectory, or in a 
different geometry (scan range, or distance to detector, etc.) than in the 
exemplary system shown, or multiple tubes and/or detectors may be employed. 
[0019] In one embodiment of the present invention, the reconstruction 
algorithm comprises a number of steps. First, for each projection angle, the 
corresponding new image is formed as a linear combination of appropriately 
filtered versions of all, or a subset of all, acquired projection images. In this 
regard, the filters used will generally be different for each projection image to 
be filtered, and each projection angle. Then, a reconstruction of the imaged 
volume is completed from these new images by using simple backprojection. 
Each of the new images created in the first step is backprojected according to 
the associated projection angle/focal spot position. 
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[0020] If the acquisition geometry is (at least approximately) known 
beforehand, the filters that are used can be computed in advance. In addition, 
the filters can be chosen such that the obtained 3D reconstruction is essentially 
equivalent to a reconstruction obtained with a predetermined or pre-specified 
number of iterations in an equivalent standard ART reconstruction (or, such that 
the obtained 3D reconstruction is equivalent to the result obtained with standard 
ART after convergence). The filters can also be chosen to approximate such a 
solution, while having additional desirable properties (e.g., being 
computationally even more efficient, or offering improved noise suppression). 

[0021] The present invention has a number of significant benefits, and is 
extremely flexible offering a wide variety of possible modifications. With the 
present invention, only a single backprojection operation per projection angle is 
required, without any re-projection. With ART, typically a large number of 
iterations are necessary in order to achieve a high degree of accuracy. With the 
present invention, however, only a modification of the used filters is required 
(according to the number of iterations used in ART) to obtain a reconstruction 
that is equivalent to a result obtained with a predetermined number of iterations 
used in ART. Generally, to obtain satisfactory results with the present 
invention, the filters will correspond to a relatively high number of iterations in 
standard ART. Thus, the present invention reduces the computational effort and 
can increase the speed of the computation significantly. 

[0022] The present method is also robust with respect to inconsistent 
projection images. In standard ART, inconsistencies between different 
projection images can lead to "non-convergence" (e.g., limit-cycle type 
behavior) of the iterative reconstruction. Because the present inventive method 
is non-iterative, this problem does not arise. 

[0023] The "exact version" of the approach according to the present 
invention requires N filtering steps for each of the N projection images (where 
N is the total number of acquired projection images). In addition, with the 
present invention, an approximation of the "exact" method can be obtained that 
speeds up the computation by replacing a filtered version of one projection 
image by a filtered version of another projection image, with a correspondingly 
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modified filter. By doing this, the number of filtering operations can be reduced, 
since different filters acting on the same image can be combined into one. In 
general, the modification of the filters will consist of a simple shift/translation 
such that the position of the filter with respect to the projected central slice of 
the volume to be reconstructed remains constant. This is depicted in Figure 2 
where a filtered version of a projection image acquired with detector 22 and 
corresponding to focal spot A can be approximated by a filtered version of the 
projection image corresponding to focal spot B. In Figure 2, the volume to be 
reconstructed is identified by reference numeral 26 and the central slice of the 
reconstructed volume is identified by reference numeral 28. By taking an 
arbitrary point 25 located on central slice 28, and by drawing logical lines 23 
and 24 through the respective focal spot locations A and B and point 25, one 
can identify the distance A by which the projection of the point is shifted in the 
respective projections. Accordingly, a modified filter can be created by shifting 
the filter by distance A. Filtering the projection image corresponding to focal 
spot B with the modified filter can serve as an approximation of filtering the 
projection image corresponding to focal spot A with the original filter. In the 
limit, this type of approach can be utilized to design filters for a filtered 
backprojection-type reconstruction in tomosynthesis, where every projection 
image is filtered (only once) and then backprojected, to form the reconstructed 
dataset. It is to be recognized that an alternative shift A can be determined with 
respect to a region of interest within the imaged object that does not coincide 
with the central slice 28 It is to be recognized, that there are also other ways to 
modify a filter such that the filtered version of one projection image can be 
approximated by a filtered version of another projection image with the 
modified filter. 

[0024] The present reconstruction method can secure improved image 

quality in the reconstruction since there are no intermediate steps / iterations as 

in ART. The iteration steps in standard ART require interpolation of the image 

data and the data undergoes the associated potential loss of image detail in the 

backprojection and re-projection steps at every iteration. ART iterations also 

include a discretization in the z component/height, thereby potentially 
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introducing artifacts in the reconstruction. The mechanism that leads to these 
artifacts is circumvented in the algorithm according to the present invention. 
[0025] With the present invention, in addition to the original projection 
images, only the new images (i.e., the combinations of filtered versions of the 
projection images) need to be stored. This reduces the data storage 
requirements. Furthermore, if only a part of the full volume needs to be 
reconstructed, this can be accommodated. In ART, in contrast, in each iteration 
the whole volume needs to be reconstructed, stored, and re- projected, for use in 
the subsequent iteration step. 

[0026] Since only a single backprojection step is utilized in the present 
invention, it can be combined with, for example, order statistics-based 
backprojection (OSBP) techniques for minimizing artifacts. OSBP could 
replace the simple backprojection step, but will in turn increase the 
computational burden in the backprojection. This makes it impractical for use in 
combination with standard ART, where many backprojection steps are required. 

[ 0027 ] A reconstruction obtained with the "exact version" of the present 
method (which means that none of the "approximations" is used, which can 
speed up the computation), satisfies the re-projection consistency constraint, 
i.e., an iteration does not lead to a further improvement. However, when used in 
combination with OSBP, or some types of constraints (e.g., thresholding at pre- 
defined minimum and maximum values in the reconstruction) an iterative 
update using the present method can be applied with extremely few iterations to 
achieve convergence. In one embodiment, the constraint is applied after a first 
reconstruction is performed. Then, the reconstructed dataset is re-projected, the 
difference to the original projection images is taken, reconstructed, and then 
used to update the previously reconstructed dataset. This update can then be 
iteratively repeated. 

[0028] The present reconstruction method offers the flexibility to 
partition the projection angles/focal spot positions into different subsets, to 
perform a reconstruction (according to the present invention) with respect to the 
projections associated with one of the subsets, and then iterate by alternatingly 
going through the different subsets, and updating the reconstructed volume. If 
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the subsets of focal spots are sets of one, then the resulting iterative procedure is 
exactly the standard ART reconstruction strategy. If only one "subset" exists 
(which contains all the focal spots), then a one-step solution can be obtained, as 
set forth above. In other instances, "intermediate" cases can be obtained and 
used, for example, to trade off the computational load of each iterative update 
with the number of iterations needed, which makes it more efficient to exploit 
additional constraints. 

[0029] Figure 7 is a flow diagram which illustrates another embodiment 
of the present invention. This embodiment is referenced generally by numeral 
100 and includes an iterative update with a partition of the projection images 
into one or more sets of projection images. 

[0030] Ina first step, the projection images 110 are partitioned into one 
or more sets of projection images 120. Each projection image should be 
contained in at least one such set of projection images. Then one of these sets 
"A" is chosen 130, and a reconstruction is performed based on the images 
contained in this set 140. It is to be recognized that in the case where the 
projection images are partitioned into only one set of projection images 
containing all the images, then the reconstruction step is a one-step 
reconstruction as described above. An (optional) constraint is then applied 150 
to the reconstructed three-dimensional information 160. The three-dimensional 
information is re-projected 170 and the difference images are then computed to 
the true projection images 180. A different set of images can now be selected 
and a reconstruction based on the corresponding difference images can be 
preformed 190. The thereby reconstructed three-dimensional information is 
used to update the three-dimensional information 200 obtained in the previous 
iteration step, to yield an improved three-dimensional dataset. By repeating 
these steps, further improved three-dimensional datasets are obtained. 
[0031] Note that in these steps the different sets of projection images are 
typically cycled through so that all projection images contribute repeatedly to 
the iterative update step. It is also to be recognized that, if the sets of projection 
images are sets containing one projection image each, the iterative 
reconstruction as set forth corresponds to a standard ART reconstruction 
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strategy. The iterative update using a partition of the projection images into sets 
of images is particularly useful in situations where a constraint is applied to the 
reconstructed dataset, followed by an iterative update. In this scenario, the 
computational load of each reconstruction step is smaller (as compared to the 
iteration using all projections), while the overall image quality may be 
comparable. 

[0032] Since the filters used in the reconstruction according to the 
present invention are directly accessible (unlike in ART, where they are only 
implicitly defined), the filters can be modified, for example, in order to achieve 
a trade-off between reprojection consistency, and noise propagation. This may 
be particularly useful in cases where the amount of noise differs between 
projection images (e.g., due to a different x-ray technique), where for example 
only the filters applied to the "noisy" images need to be modified. Also, only 
the high-frequency characteristics of the filters can be utilized to obtain a 
reconstruction of the fine details (i.e., high frequency content) in the imaged 
volume. In the situation where the present method is combined with OSBP, the 
filters can be modified to counteract some of the expected "loss in reprojection 
consistency" due to the OSBP operations. (With OSBP, backprojected values 
are generally not taken into account at every height in the reconstructed volume. 
Larger structures, corresponding to lower frequencies, are typically taken into 
account at a larger fraction of the height of the reconstructed volume than 
smaller structures, corresponding to higher frequencies. The filters can be 
adjusted to counteract these effects to some degree, thereby maintaining a 
"better" re-projection consistency). 

[0033] For a linear tube trajectory, for example parallel to the columns 
of the detector, the filters are one-dimensional (ID), acting on columns of the 
image, and are identical for all columns of the image. Generally, in cases where 
the tube trajectory is linear with a constant height above the detector plane, the 
filters that are derived according to the present invention are 1-D filters, acting 
on lines of the images that are parallel to the tube trajectory. 
[0034] In situations where the x-ray tube trajectory is essentially linear, 
the filters according to the present invention can be combined with an additional 
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transversal filtering component, in order to further enhance the reconstructed 
image quality. 

[0035] Also, the reconstruction method of the present invention can be 
used for quantitative reconstruction, since it satisfies the re-projection 
constraint. 

[0036] The filters that are used in the reconstruction method of the 
present invention can be derived in various ways. In particular, in one 
embodiment, the filter construction will proceed in iterative steps, where each 
step mimics an iterative update step in a corresponding ART reconstruction. 
The filter derivation is motivated by the fact that a simple backprojection (of a 
single projection image) into a volume of constant thickness, followed by re- 
projection (with respect to a different focal spot position) generates an image 
that can also be (approximately) generated from the initial projection image by a 
linear filtering step with an appropriate one-dimensional (ID) boxcar filter. The 
combination of a backprojection step followed by a re-projection step (acting on 
a single projection image) can therefore be replaced by a simple convolution 
with a suitable filter/kernel. 

[0037] Additive ART essentially consists of a sequence of iterative 
backprojection/re-projection steps, combined with appropriate add/subtract 
operations. With the present invention, all these operations are collapsed into 
the construction of a set of appropriate filters such that filtering the projection 
images and suitably combining the filtered images, followed by a simple 
backprojection, is equivalent to additive ART. For every projection image a 
new image is formed as a (weighted) sum of appropriately filtered versions of 
all acquired projection images. Then, a backprojection step is performed using 
these new images, in order to obtain a reconstruction that is equivalent to one 
obtained with additive ART. This approach is also extremely flexible. For 
example, computationally less expensive approaches that approximate this exact 
result can easily be developed within the present framework. 
[0038] In tomosynthesis, the focal spot positions are typically assumed 
to be (at least approximately) at a constant height above the detector. From this 
assumption it follows that, for any given point located in a plane parallel to the 
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detector, a line connecting the projections of this point with respect to two 
different focal spot positions onto the detector has a predetermined length and 
orientation, both depending only on the height of the plane, and the relative 
location of the considered focal spots, but not depending on the specific point 
within the plane. This is shown in Figure 3 and identified generally by the 
reference numeral 30. Using now the fact that in ART reconstruction for 
tomosynthesis, the volume to be reconstructed will typically be limited by a 
minimum and a maximum height above the detector, and the fact that the 
projection of a line is a line, it follows that backprojection with respect to one 
focal spot, followed by re-projection with respect to another focal spot, is 
(approximately) equivalent to a convolution with an appropriate ID boxcar 
filter. 

[0039] As shown in Figure 3, for any point at a given height above the 
detector 34, the projections with respect to two focal spots A and B are at a 
constant distance from each other, independent of the specific point chosen. For 
example, the projections of any point at height hi are at a distance of Al from 
each other. It follows that a backprojection with respect to focal spot A, into the 
volume 36 between heights hi and h2, followed by a re-projection with respect 
to focal spot B, is (approximately) equivalent to a convolution with a boxcar 
filter. The boxcar filter has a length of (A2-A1), and an offset of (Al+A2)/2 (as 
compared to a "centered" boxcar filter of the same length). 

[0040] As indicated above, the construction of the filters relies on the 
fact that, together with an assumed minimum and maximum height above the 
detector for the volume containing the imaged object, a backprojection of a 
single image (associated with one projection angle), followed by a re-projection 
(with respect to a second projection angle) is equivalent to filtering the image 
with a boxcar filter of a given length and offset, and with integral one. The 
exact filter, however, is not a boxcar filter, but has a non-constant amplitude. 
However, a boxcar filter can generally be used as a good approximation. 

[0041] A standard ART strategy is, for example, to perform a simple 

backprojection based on a single projection image. The resulting 3D dataset 

represents the initial reconstruction. The reconstructed volume/dataset is then 
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re-projected with respect to a projection angle corresponding to a second 
projection, and the difference between the true second projection image, and the 
re-projected image is computed. This difference image is then backprojected 
(according to the projection angle corresponding to the second projection 
image) and added to the initial reconstructed volume. The quality of the 
reconstruction is improved by iteratively using the same steps: re-project the 
current reconstruction (with respect to some projection angle), backproject the 
difference between true projection image (for that projection angle) and re- 
projected image, and add to the current reconstruction to obtain an updated 
reconstruction. Formally, this strategy can be written in the following form, 

0<"> =0^ + B i(n) (l Hn) -P Hn) (O in -»)) (1) 

where 0 (n) denotes the estimated reconstructed object at the n-th iteration, Ij (n ) 
denotes the projection image corresponding to projection angle i(n), and Pi (n) 
and Bi(n) are the associated projection and back-projection operators, 
respectively. The initial condition is usually chosen to be O (0) = 0. 
[0042] There are various ART strategies, which differ mainly in the 
order in which the different projection images are visited, the chosen weight for 
the update term, and whether each update is based on a single, or multiple 
images (in the multiple image case, the backprojected images, or backprojected 
differences, are generally averaged). However, these different strategies are 
aimed solely at improving the speed of convergence, and they all converge to 
the same resulting reconstruction. For this reason, any one of these strategies 
can be used to derive the filters according to the present invention. 
[0043] To explicity derive the filters, the specific strategy of the 
example given above can be followed, in which case one obtains (by rewriting 
Cr " according to the iteration Equation (1)) an iterative update of the form: 

0<«>=0<«- 2 > + b h ^ } (V ir ^,,-(0 M )) + 

B i(n) (Au) - P i(n)(° (n ~ 2) +^(„-I,( / /(„-l) - P i(n-\)(° (n ~ 2) ))) 

[0044] As stated above, backprojection (with respect to one projection 
angle), followed by re-projection (with respect to a second projection angle) is 
essentially equivalent to filtering with the corresponding boxcar filter. If the 
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filtering operator corresponding to a backprojection with respect to angle i, 
followed by a reprojection with respect to angle j is denoted with K(ij), then the 
previous equation can be rewritten as: 

O"" = 0 ( "- 2 > + £*„„_„ ( /„„_„ - P i{n _ ]) 0 { "- 2> )) + 

*«■> ((/«., " ^ n ,(0 ,n - 2) + B i{n _ n (I iin _ u - /^ ( „_ 1) (C• ( "- 2, ))) = (3) 
C ( "- 2, + 5 i(n _ 1) (/, n . 1) -^ (n _ 1) (0"'- 2) )) + 

*«■>('«., -i? ( .,O ( " 2, -«0X«)wX«-l))/«.-, ) + ^(iX»)wXii-i))Vii(^"" 2> »^ 

[0045] By successively expressing the reconstructions 0 (n 2) etc. by 
''previous" reconstructions, and by rewriting any combination of back- 
projection/re-projection PjBi as the corresponding filtering operator K(i,j), 
together with the initial condition O (0) = 0, the final equation is obtained where 
0 (n) is expressed as a linear combination of backprojected images, each image 
consisting of a linear combination of appropriately filtered versions of the 
original projection images. The involved filters are expressed as the (weighted) 
sum of combinations of the elementary filters K(i, j). 

[0046] By letting the "iteration number" n increase (for example to 
n=10N where N is the number of projection images - a large n ensures 
convergence of the equivalent ART reconstruction), the filters F can be obtained 
that are required to compute the reconstructed dataset: 

&*=I. i B i [L J F w {iJ)I J ) (4) 

which is essentially identical to the result obtained after n iterations of additive 
ART (using the identical update strategy as outlined above). 
[0047] The above suggests to create the filters F by successively 
substituting 0 (n 2) etc. in Equation (1), thus obtaining a polynomial in K(i,j) for 
each of the filters F, and then explicitly computing each F from the "elementary 
filters" K(ij). Using this approach, and using the approximation of the filters K 
as boxcar filters, the filters F can also be explicitly expressed as piecewise 
polynomials, or as B-splines. 

[0048] A second strategy to derive the filters F, for example, can be to 
assume a representation of Equation (4) for 0 (n) , and to use the update Equation 
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(1), which will provide a recursive method to compute the filters F(i,j). 

Specifically, we have: 

0M =0 <*-.> +Bi(a) (l iln) -P i(n) 0^)) 

= Z,B, ) + £,,„• ( / i(n) - P^ftB, (2 ; F"-"( f -, )) 

= (S .F""'>(i, ,)/. ) + fi i(n) (/, n) - Z^ (n) 5,. (Z y F<-"(/, ,)/, )) 
= Z,S,.(Z.F<"-%^ (5) 

and it follows that: 

F ( ->(i, j) = F ( -"(i f 7) /or f * /<«), (6) 

F (n) (/, y) = F ( -" !) (i f j)-Z k K(i,k)F {n ~ l) (k, j) for i = i(n), j * i(n), (7) 

F {n \iJ) = F (n - l \i, j) + S-X k K(i,k)F {n - l) (kJ) for i - i(n)J = i(n). (8) 

where 5 denotes the Dirac impulse. 

[0049] As noted above, different ART update strategies can be used as 
well, leading to modifications of these recursive definitions of the filters F. 
Also, since all additive ART strategies converge to the same end result (i.e., the 
minimum norm solution to the projection equation), there exists a set of filters 
that correspond to that solution, and it is possible to derive these filters directly, 
for example by using a fixed point condition for the recursive update equation 
for the filters. This condition leverages the fact that, once the iteration has 
converged, an additional update step will leave the filters unchanged. 

[0050] One way to derive the filters F(i,j) which can be used in the 
method of the present invention to compute a reconstructed volume that is 
equivalent to a reconstruction obtained with ART after convergence is to utilize 
the re-projection consistency constraint that the corresponding reconstruction 
(theoretically) satisfies. This constraint says that the reconstructed volume, 
when reprojected with respect to one of the considered focal spot positions, 
results in a projection image that is identical to the true projection image for that 
focal spot position. Using a presentation as in equation (4), this constraint can 

be expressed in the following form, I k = P k O = P k W^iW^ JV s )) (f° r a ^ 
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k). Writing this expression (for all k) now in matrix form, one obtains 



(i 


r \ ( 
i 'i 

• # 

• • 

p 

* 
• 


■fa - B N y 

••• K(l,N)^ 

m 

• 


' f(U) •• 

• ■ 

• • 
• 

( F(l,l) • 

• « 
• 


• F(l, AO ^ 

• 


• 

VA/> 
♦ 






m 

K(N,N)j 


• 

,F(AU) 


• | 

F(AT, AO, 


• 



and it follows that the filters F(i j) can be determined as the entries of the matrix 
that is the "inverse" of the matrix containing the elementary filters K(i,j). One 
way to obtain the filters F(i,j) uses a Fourier space representation of the filters 
K(i,j) and F(i,j). Using this representation, one obtains a separate equation 
(corresponding to the previous equation) for each frequency in Fourier space, 
and the filters F(i,j) (in their Fourier space representation) can be easily 
determined by computing for each frequency the inverse (or, generalized 
inverse, or Moore-Penrose inverse, for example with the help of a singular value 
decomposition) of the matrix containing the corresponding Fourier coefficients 
of the filters K(i,j). It is to be recognized, that other approaches to determine the 
filters F(i,j) that will lead to a reconstruction that is equivalent to a 
reconstruction obtained with ART after convergence can be used as well. 
[0051] Figure 4 depicts examples of filters F(i,j), for a tomosynthesis 
scenario where five different projections are acquired and used for the 
reconstruction, with a tube trajectory on a straight line, parallel to the columns 
of the detector. In this scenario, all filters are ID, and aligned with the pixel 
columns of the projection images. In this scenario, a new image "i" will be 
computed as the sum of filtered versions of all projection images "j", where 
each image "j" is filtered with the corresponding filter F(i,j). The reconstructed 
volume is then obtained, for example, by backprojecting and averaging these 

new images ^ F(i, j)I . . It should be appreciated that, while the elementary 

filters K(i,j) are ID filters, the filters F(i,j) will typically be 2D (two- 
dimensional) filters. Only in the case when the tube trajectory is linear (i.e., the 
tube positions for all projections are located on a line) and at a constant height 

above the detector plane will the filters F(i,j) be ID filters. 
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[0052] As mentioned above, from the filters F(ij), filters can be derived 
that can be used in a more standard filtered backprojection type reconstruction, 
where each image is only filtered once, and then backprojected. One example 
of such filters derived from the filters shown in Figure 4, is shown in Figure 5. 
The filters shown in Figure 5 correspond to a standard filtered backprojection 
type reconstruction and were derived from the filters depicted in Figure 4. In 
this scenario, only a single filter corresponding to each projection image is 
depicted. After the filtering step, the volume of interest is reconstructed using a 
backprojection operation. 

[0053] As mentioned above, suitable boxcar filters represent generally 
satisfactory approximations of the elementary filters K(i,j). If even better 
accuracy of the reconstruction is required, the exact filters K(i,j) can be used. 
Figure 6 illustrates the derivation of an exact filter corresponding to a 
backprojection with respect to one focal spot, followed by a re-projection with 
respect to a second focal spot. Simple backprojection of a single pixel (the 
pixel is indicated by the interval of length 5) with respect to the focal spot 
position A assigns to each location within the intersection of the 
"backprojection cone" and the volume 52 to be reconstructed a constant value 
(given by the image value at that pixel, divided by the height/thickness of the 
volume of interest). It is assumed that the image values are normalized with 
respect to the pathlength of each ray through the volume of interest. That is, in 
the normalized image domain, a ray through a homogeneous material of 
constant thickness will generate a constant image value, independent of the 
location of the pixel, or the incident angle of the specific ray. With these 
assumptions, the re-projection (with respect to focal spot B) of this 
backprojected pixel, at a given location which is defined by the offset "x" from 
the original pixel is the difference h2-hl, i.e., the height difference between the 
points where the ray enters and exits the backprojected "cone" associated with 
that pixel. 

[0054] From Figure 6, the specific elementary filter K can be 
determined. Specifically, it can be seen that 
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*2 = 



(S + x)H 



, and/z, = 



xH 



, and therefore 



A + J + jc 



A + jc 



//(£ + *) (A + ;c)-//jc(A + c? + jc) 
(A + £+ jc)(A + ;c) 

//<?A 
( A + S + jc)( A + x) 

//<?A 
(A + x) 2 + <?(A + ;t) 



where the last line is an approximation which is valid in the case that 5 is small 
as compared to A, which is generally true. This previous equation defines the 
kernel K exactly, for all "x" such that hi and h2 are inside the volume of 
interest. It is also clear from Figure 5 how the "filter boundaries" can be 
derived. Obviously the kernel K is zero for all 'V such that neither hi nor h2 
are contained in the volume of interest, while both "ramp-up" and "ramp down" 
of the filter are exactly one pixel wide. Furthermore one can see that, in the 
case where the volume of interest is "close to the detector", i.e., for "x" small as 
compared to A, the boxcar filter is indeed a good approximation of the filter K. 
It is also obvious from this derivation that the filter K is shift-invariant, which is 
then also true for the filters F. 

[0055] While the invention has been described in connection with one or 
more embodiments, it is to be understood that the specific processes and 
procedures which have been described are merely illustrative of the principles 
of the invention, numerous modifications may be made to the subject matter 
described without departing from the spirit and scope of the invention as 
defined by the appended claims. 



(7) 



(A + jc) 2 
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